How efficient is towing a cargo by a micro-swimmer? 
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We study the properties of arbitrary micro-swimmers towing a passive load through a viscous 
liquid. The simple close- form expression for the dragging efficiency of a general micro-swimmer 
dragging a distant load is found, and the leading order approximation for finite mutual separation 
is derived. We show that, while swimmer can be arbitrarily efficient, dragging efficiency is always 
bounded from above. It is also demonstrated, that opposite to Purcell's assumption, the hydrody- 
namic coupling can "help" the swimmer to drag the load. We support our conclusions by rigorous 
numerical calculations for the "necklace-shaped" swimmer, towing a spherical cargo positioned at a 
finite distance. 

PACS numbers: 



In the recent years there has been an increasing inter- 
est in propulsion at low Reynolds numbers, both theo- 
retically [r, ^, BM, i, 0, S, [3, E [13, M, Hi and experi- 
mentally [6, 8,ll6|. These and other works have improved 
our understanding of the basic properties of locomotion 
on small scales. However, it is not sufficient to under- 
stand the mechanisms and properties of free microswim- 
mers alone - it is necessary to estimate the performance 
of these swimmers as propellers that tow a useful cargo, 
e.g. therapeutic load or miniature camera. This ques- 
tion, which attracted only limited attention, had already 
shown to have some non-trivial answers: E. M. Purcell 
had studied the particular case of a rotating helix 
pushing a spherical particle under the assumption of neg- 
ligible hydrodynamic interactions. He showed that due to 
the structure of grand-resistance matrix, which connects 
the force and torque on a body to its translational and 
angular velocities, the optimal rotating-propcUer should 
have the same size as the load. 



In this paper we address arbitrary shaped swimmers 
and cargoes, and investigate the effect of their mutual hy- 
drodynamic interaction on performance of the swimmer 
as load propeller. We find that, while propellers that can 
enclose a load within may theoretically have arbitrarily 
high efficiency (consider, for instance, the "treadmiller" 
) , the dragging efficiency of a swimmer towing a remote 
load is always bounded from above, and there is an op- 
timal propeller/load size ratio. We also find, in contrast 
to Purcell's assumption that in many cases the hy- 
drodynamic coupling between the load and the propeller 
enhances the dragging efficiency. We support our theory 
by numerical calculations for a necklace-shaped swimmer 
ll\ towing a spherical cargo. 



We start our analysis by considering a micro-swimmer 
(i.e. a propeller) dragging a distant load. In this case, 
we can neglect the mutual hydrodynamic interaction and 



calculate the dragging efficiency as 

KiVi 



Pd 



(1) 



where Ki is the resistance coefficient of the load [21| , Vd 
is the dragging velocity and Pd is the power expended by 
the swimmer to drag the load with velocity Vd- We also 
define the propeller's efficiency in the same fashion, Ss — 

K 

p " , where Ks is the swimmer's resistance coefficient, 
Vs is the speed of the unloaded propeller (at the point 
where the load is anchored) and Ps is the power expended 
in swimming without load. 

Note that in a general case of swimmer propelled by a 
sequence of geometrically non-reciprocal periodic strokes 
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(e.g. three- link Purcell's swimmer [s Mla . Il7l |. push- 
mepuUyou \4, three-sphere swimmer shape defor- 
mations p. Il9| and others), the swimming efficiency is 
conventionally defined using stroke-averaged quantities 
[l^ . However, since max | ^p^° | > '^"^j^^^ (where () 
stands for average over a stroke period; the maximum is 
taken over the stroke period), the maximum of ([T]) over 
a stroke period is an upper bound for the conventional 
efficiency. In the case of swimmer propelled without the 
shape change (e.g . rotating flagella 14 1, treadmiller 0], 
twirling torus and others), the two definitions coin- 
cide. They are also practically equivalent for swimmers 
performing small-amplitude strokes, with Kg w Const 
Also, note that llj is not just the swimming effi- 
ciency re- written for "swimmcr+load" as a new swimmer, 
since we aim to compare the power expended in dragging 
the load alone by swimmer and by an external force. 

We will now calculate the dragging efficiency for swim- 
mer characterized by a resistance coefficient and 
swimming efficiency e^, dragging a load characterized by 
a resistance coefficient Ki , which we will assume are both 
not rotating (it is known [l[ that a rotating swimmer 
is less efficient than a non-rotating one). By Lorentz 
reciprocity [lo|, if (vj^ajk) and (vpa'.j^) are the velocity 
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and stress fields for two solutions of the Stokes equations 



djffij = in fluid domain E then 



as 



Tij dSj — 



(2) 



Using (21) with {vi , aij ) being the velocity and the stress 
fields for a swimmer dragging a load, and {v^,a[j) being 
the velocity and stress fields for the "unloaded" swimmer 
and the load co-dragged by the external force with the 
swimmer's velocity, we readily obtain 



Pd = Ps+i - {Vs - Vd) ■ Fs+i 



(3) 



Here Pd is the power expended by the swimmer to drag 
the load, Pg+i is the power dissipated by viscosity in the 
case of "unloaded" swimmer and the load co-dragged by 
the external force, Vs is the velocity of the free swim- 
mer, Vd is the dragging velocity and Fg+i is the force 
required to tow the load with velocity Vg. Vd can be 
found by equating the sum of the viscous drag forces on 
the swimmer and the load to zero. Exploiting the lin- 
earity of Stokes equation and neglecting hydrodynamic 
interaction, we get: 



Vd 



Ks + Ki 



(4) 



As expected, Vd goes to zero for infinity large load and to 
the swimmer velocity for a vanishingly small load. Ne- 
glecting hydrodynamic interaction, we can use i^^+z — 
-Vs Ki and P^+i = Ps + V^ Ki, that together with g]) 
and dragging power reads 



Pd = P. 



Vf 



KiKs 
Ki+Ks 



(5) 



For small loads, Ki ^ Kg, Eq.® gives the power of 
the free swimmer plus the power of dragging the load, 
and for large loads, Ki :s> Kg, this give s the power of 
an anchored swimmer (i.e. a "pump" Substituting 
dS]), the swimmer efficiency and ^ into ([T]) gives 



£d 



(r + l)(r±i 



+ r 



(6) 



where r — Ki/Kg. The dependence of the dragging ef- 
ficiency Ed on Es and r is plotted in Fig. ([T|). Eq. Q 
shows that unlike the swimming efficiency, which, for 
some swimmers, can be arbitrarily high 0,3, the drag- 
ging efficiency is bounded by Ed < 7^ < 1 even for 
Es = 00. This means that enclosing a cargo within the 
swimmer can be much more efficient than towing a re- 
mote one, and that there is an optimal swimmer size 
for any swimming technique (including swimming tech- 
niques in which r is varying periodically). As one might 
expect, Ed is a growing function of Eg- However, while 
for inefficient propeller (like a rotating helix) the opti- 
mal size is about the same as the load size, the efficient 




FIG. 1: (Color online) The dragging efficiency, Sd, as function 
of the propeller's efficiency Sa and the size ratio Ki / Ks 



swimmer with Eg ^ 1 (e.g. "pushmepullyou" [J|) will 
be efficient as propeller only if it is much larger then the 
load. Thus, the naive intuition saying that the swim- 
mer's efficiency alone controls the dragging efficiency, is 
not always right: in some cases less efficient but larger 
propeller is advantageous. 

Now let us estimate the effect of hydrodynamic inter- 
action between the propeller and the passive cargo sepa- 
rated by distance d. For finite separation distances it is 
no longer valid to assume that -Fs+z = —Vs Ki. However, 
the force must still be linear in the dragging velocity and 
we can write Fd = XiKi Vd- In the same way, the force 
on the swimmer must be proportional to the changes in 
the velocity, so we will denote it by Fs — \s{Vs — Vd)Ks- 
Since the forces must still sum up to zero, the dragging 
velocity is 



Vd^ 



VKs 



Ks 



Ki 



(J) 



Comparing the velocity in ([7]) to that for infinite separa- 
tion it can be readily seen that the deviation between 
the two depends on the ratio ^: if > 1 the velocity 
will be lower then that in and if < 1 the veloc- 
ity will be higher than the infinite distance case. Since 
generally A^, A; < 1 and for asymmetric configura- 
tions the resistance coefficient of the larger object will be 
almost constant, we can conclude that a large swimmer 
will drag faster when positioned close to the load, while a 
small swimmer will drag faster when located far from the 
load. 

Assuming that the separation distance is large enough, 
so d > max{i?/, where Ri and Rs are the hydrody- 
namic radii of the load and the propeller, respectively, 
we can now estimate the power needed for the swimmer 
to drag the cargo: it is known 15| that for any swimmer 
Ps — Pp — Pg, where Pp is the power needed by the pump, 
i.e. the anchored swimmer, Pg is the power needed by the 
swimmer when it is swimming freely and Pg is the power 
needed to drag a "frozen" (immobile) swimmer with the 
swimming velocity. If we use this relation by treating the 
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swimmer plus the load as a modified swimmer, we can 
estimate the power needed to drag the load. In this case, 
Pg is just the power needed to drag both the load and 
the (frozen) swimmer with velocity Vd provided by ([7]), 

where Pg(s) = V^K^ is the 



which is Pg = Pg{s)^rE. 



+1 



power needed to drag the immobile swimmer. As an ap- 
proximation, we will assume that the power expanded by 
the pump does not depend on the proximity of the load, 
since Pp = Pp{s) + ^'[("f )^] (^p(s) is the power expanded 
by the pump when the load is absent). Together, this 



gives 



Pd^Ps+ V.K, 1 



A. 



(8) 



where Pg is the power expended by the swimmer when 
the load is absent. Obviously, Pd > Ps and the equality 
holds only when Ki — 0. 

Substitution of (O and in H]) yields 



£d 



A, r 
As 



1 + ^ -A 



Aj r 
As 



(9) 



where r — Ki/Kg. For As = A/ = 1 (H]) reduces to (O, 
as anticipated. Comparing the efficiency in ^ to that in 
^ one can conclude that in cases where > 1 (i.e. the 
swimmer is smaller than the load) , the efficiency is lower 
when the hydrodynamic coupling is not negligible, and it 
would better be separated from the load. If the swimmer 
is much bigger than the load, which implies As — 1 and 
< 1, the efficiency is higher than in the case with 
no coupling! Thus propeller bigger than the load should 
be positioned closer to the load, opposite to Purcell's 
assumption [3]. Eq. ^ also tells that the efhciency is 
bounded by which for large propeller towing a small 
load can theoretically be grater than 1. However, we 
could not find such an example. 

We can now estimate As and A; as functions of d, 
using the Oseen tensor [lO|. As the first order ap- 
proximation, we will assume both the swimmer and the 
load can be modeled as spheres with hydrodynamic radii 
i?s ~ and Ri = respectively [2^. In this case, 
it can be readily shown that for d ^ max{i?;,i?s} [lot : 

2{2d^-3dRi) _ 2(2d^-3rffi3 
A' 



A 



4d--9R,R. ' = 4d-^-9R:Rj - Substituting these 
expressions into ([7]) gives the expression for the dragging 
speed, 



Vd = Vs 



26-37 



2{S-3r + 5r) ' 



(10) 



where 5 = d/ Rg is the scaled separation distance. Simi- 
larly, substitution of A; and As into ([9]) yields 



£d 



2S-3 , 
25-3r ' 
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2(2(5^-3ri5) 



4i52-9r 



25-3 „ 
2(5-3r' 



(11) 




FIG. 2: (Color online) Schematic of the necklace- like propeller 
towing a spherical load. The arrows show the direction of the 
rotation of spheres in the propeller; the propeller is pushing 
the load in front of it. 



Expanding ([TT|) for small | gives £d ~ ^d(oo) + 
[l — (1 -I- £s)r^] -I- . . . , where £^(00) corresponds to 



the no-hydrodynamic-interaction approximation for the 
dragging efficiency ^ . The 1 /(5-term in the above expan- 
sion shows that for r > ^^]^^ the dragging is retarded 
in comparison to the infinite separation result ([6]), i.e. 
£d < £d{oo), while for r < ^7^=^, the dragging is en- 
hanced due to the hydrodynamic coupling, as Ed > £d(oo)- 
Interestingly, r — ^^^^^ corresponds to the maximum of 
£d(oo)- However, it is not the optimum of e^, which shifts 
to higher values at smaller r's. 

We shall now test the proposed theory for the load 
dragged by a rotary propeller. Imagine the necklace- like 
ring (see Fig[21) of iVp = 8 nearly touching rigid spheres 
(separated by the distance of 0.05a) of radius a. The 
necklace lies in the xy plane and in a cylindrical polar co- 
ordinate system (z, r, ip), each sphere rotates at the con- 
stant angular velocity H — we^, which, in the absence 
of external forces, causes the necklace to swim along the 
normal to the plane of the necklace in the positive z direc- 
tion |lll | . Performance of this swimmer as cargo propeller 
is tested for a spherical particle positioned at arbitrary 
distance along z-axis [23|. The distance that separates 
the plane of the propeller (z — 0) and the load's surface 
is denoted by . 

We use the Multipole Expansion method [ll| and 
construct the rigorous solution of the Stokes equations 
as su perp osition of Lamb's spherical harmonic expan- 
sions [10|. The no-slip conditions at the surface of all 
spheres are enforced via the direct transformation be- 
tween solid spherical harmonics centered at origins of dif- 
ferent spheres. The accuracy of calculations is controlled 
by the number of spherical harmonics, L, retained in the 
series. The truncation level of L < 7 was found to be suf- 
ficient for all configuration to achieve an accuracy of less 
than 1%. The dragging efficiency ([T]) for this particular 
swimmer reads 



ed 



NpTu 



(12) 
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FIG. 3: (Color online) Numerical results for the "necklace- 
shaped" propeller made of 8 co-rotating spheres of radius a 
{Vs/auj ~ 0.316, Es — 0.0339, Rs — 3.083a), towing a spheri- 
cal load of variable size located at: d, — (magenta), d* = a 
(blue), d, = 4a (red) and d* = 10a (yellow); the solid line 
corresponds to the infinite-separation result Q; the dashed 
lines are the far-field approximations for d* = 10a. (a) the 
scaled dragging speed, Vd/aw; (b) the dragging efficiency, Ed- 



where T = x (cr-n) dS is a hydrodynamic torque 

exerted on each (ith) sphere of the propeller towing the 
load with being the radius vector with origin at the 
ith particle center. The values of if/, T and Vd are de- 
termined numerically and the resulting scaled dragging 
speed Vd/auj, and efficiency, Ed, are plotted vs. the size 
ratio r in Figs, and (6), respectively. The agreement 
with the asymptotic results (|10lll[) (via S = d^, / Rs + r) 
is excellent for small loads (r < 1) even at moderate dis- 
tance of = 10a (i.e. S ~ 3.24). It can be readily seen 
that there is an optimal swimmer-load size ratio in all 
cases. Interestingly, Fig. shows that while for mod- 
erate separation the dragging velocity decays with the 
increase in the load size, at close proximity it may actu- 
ally be higher than the velocity of the unloaded swimmer. 
This is a direct consequence of Eq. ([7]), which does not 
assume large separation: the fluid velocity in the cen- 
ter of the "necklace" is larger than the swimming speed. 
This means, that in order to pull a load positioned at 
ds, = with the swimmer's speed, the applied force must 
act in the direction opposite to that of the velocity, so 
that A; < and Vd> Vg. 

The numerical results confirm the qualitative depen- 
dencies arising from the far-field theory: there is a criti- 
cal size-ratio r^r (weakly dependent on d) such, that for 



r < the dragging efficiency is higher than the corre- 
sponding ed(oc) and for r > Vcr the efficiency is lower than 
ed(oo); at moderate separations Vcr ^v+e ^ expected 
from the far-field analysis. The discrepancy between the 
asymptotic and the numerical results is only observed at 
r > 1, where the assumption - ^ 1 is no longer valid. 
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We treat Ki as a scalar, since the resistance matrix is 
symmetric p^ . and it is optimal to drag the load when 
the velocity is aligned in the direction of the minimal 
eigenvalue 

If either load or propeller deviates considerably from the 
spherical shape, it can be taken care of in the far-field 
resistance tensor. Here, for simplicity, we only refer to 
the far-field interaction between two spheres 
For computational simplicity, the potential hydrody- 
namic disturbance caused by the links required to con- 
nect the propeller and the load is neglected. 



